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Context. Submillimeter galaxies are a population of dusty star-forming galaxies at high redshift. Measuring their properties will help 
relate them to other types of galaxies, both at high and low redshift. This is needed in order to understand the formation and evolution 
of galaxies. 

Aims. The aim is use gravitational lensing by galaxy clusters to probe the faint and abundant submillimeter galaxy population down 
to a lower flux density level than what can be achieved in blank-field observations. 

Methods. We use the LABOCA bolometer camera on the APEX telescope to observe five cluster of galaxies at a wavelength of 870 
/jm. The final maps have an angular resolution of 27.5"and a point source noise level of 1.2-2.2 mjy. We model the mass distribution 
in the clusters as superpositions of spherical NFW halos and derive magnification maps that we use to calculate intrinsic flux densities 
as well as area-weighted number counts. We also use the positions of Spitzer MIPS 24 /jm sources in four of the fields for a stacking 
analysis. 

Results. We detected 37 submm sources, out of which 14 have not been previously reported. One source has a sub-mjy intrinsic 
flux density. The derived number counts are consistent with previous results, after correction for gravitational magnification and 
completeness levels. The stacking analysis reveals an intrinsic 870 yum signal of 390 ± 27 ytijy at 14.5cr significance. We study the 
SiA/an - Sgioum relation by stacking on subsamples of the 24 /jm sources and find a linear relation at S24 l ,m < 300 /jJy, followed by 
a flattening at higher 24 yum flux densities. The signal from the significantly detected sources in the maps accounts for 13% of the 
Extragalactic Background Light discovered by COBE, and the stacked signal accounts for 1 1%. 

Key words. Galaxies: clusters: individual: Abell 2163, IE 0657-56 , AC 114, Abell 2744, MS 1054-03 - Submillimeter: galaxies - 
Infrared: galaxies - Cosmology: observations - Gravitational lensing 
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Submillimeter galaxies form a population of high-redshift, dusty 
star-forming galaxies that are highly obscured in the visible 
and in the near-infrared, and have a spectral energy distribu- 
tion (SED) that peaks in t h e sub millimeter (submm) wave- 
band (see e.g. iBlain et aD (I2002I) for a review). Most re- 
cent searches for submm galaxies have been based on sur- 
veys of blank sky with no known large-scale structure along 
the line-of-sight. These surveys have exploited large-format 
sensitive bolometer arrays ( e.g. SCUBA on the James Clerk 
Maxwell Telescope (JCMT), ICoppin et al]l2006t LABOCA on 



X 



the Atacama Path finder Experiment (APEX). IWeifi et al.ll2009t 
AzTEC on ASTE: lAustermann et al.ll2010tlScott et al.ll2010l~the 



South Pole Telescope ( SPTl:lVieira et al.ll2010t MAMBO on th e 
IRAM 30 m telescope: iGreve et al.ll2004 iBertoldi et al.ll2007l) . 



Those maps cover large areas at a nearly uniform noise level, 
leading to a simple selection function with a constant complete- 
ness across the field. The observations showed that source num- 
ber counts increase steeply w ith decreasing flux density at mm 
and s ubmm wavelengths (e.g. lWeifi et al.ll2009t iPatanchon et al.l 
120091) . In order to probe the faint (below a few mjy) popu- 
lation of submm galaxies, several authors have taken advan- 



ters of galaxies (e.g. 


Smail et al. 


1997 


120021 Chaoma 


i et al.l 


2002bt Knudsen et al.l 


2005L 2006 


2008 


; Johansson et al. 


20 id 


Wardlow et al.l2010: Rex et al. 2009: Eaami et alj2010h. A large 



magnification, produced for example when a source lies close 
to a critical line of the lens, may make it possible to detect a 
source with an intrinsic flux density much lower than the for- 



mal root mean square of the noise of the observation. This is 
the only method of detecting such dim sources directly. Cluster 
field observations have sensitivities that vary across the map, as 
magnified sources are "lifted" above the detection limit, and the 
selection functions are therefore more complicated. The most 
comprehensive st udy to date of submm galaxies behind lensing 
clusters is that of iKnudsen et al.1 (120081) . who analyzed SCUBA 
data from 12 galaxy clusters and one blank field, resulting in an 
effective surveyed area of 7 1 .5 arcmin 2 on the sky, but an area 
in the source plane almost twice as small. Seven sources with 
sub-mjy fluxes were detected. 



The sources revealed by gravitational lensing are prime 
targets for observatio ns across the elecnomagnetic spectrum. 
ISwinbanket~aT1(l2010l) discovered a very bright submm source, 
situated at z — 2.33, with flux density 587o^ m ~106 mjy, and 
molecular line observations showed that the amount of molecu- 
lar gas is sim ilar to that in local ultr a-luminous infra-red galax- 
ies (ULIRGs, Danielson et al.l2010b . The z ~ 3.4 submm source 
studied by llkarashi et al. J 2010|) and di scovered through use of 
AzTEC at 1 . 1 mm (I Wilson et al Il2008ah has a 880 yum flux den- 
sity measured by the Submillimeter Array (SMA) of ~ 73 mjy 
and seems to be a ULIRG as well. O n the other hand, th e 
z = 2.79 galaxy behind the Bullet Cluster dGonzalez et al.l2010l) . 
with a flux density of about 48 mjy at 870 |im, is more repre- 
sentative of the normal galaxy population with an intrinsic far- 
infrared luminosity o f a few times lO 11 ^ ([Wilson et ai1l2008bl 
iGonzalez et al.l 120091: iJohansson et al.ll2010l) . Large surveys in 
the mm (SPT) and the far-infrared (Herschel) are discovering 
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bright lensed submm galaxies (IVieira et al.ll2010tlNegrello et al.l 
l2010h . 

Another way to probe the faint part of the submm galaxy 
population is to perform a stacking analysis using known po- 
sitions obtai ned from complem entary observations at another 
wavelength. iDole et al.l (120061) used the positions of sources 
detected with Spitzer Space Telescope at 24 jum to measure 
the contribution of those sources to the 70 and 160 fim far- 
infrared background, gaining up to one order of magnitude in 
depth. Greve et al. (2010) carried out a stacking analysis of the 
LABOCA submm map of the Extended Chandra Deep Field 
(ECDF) using a large sample of near-infrared detected galaxies. 

In this paper, we extend the analysis of sub mm sources be- 
hind t he Bullet Cluster recently presented in iJohansson et alJ 
j2010l) by four additional galaxy cluster fields observed with 
the LABOCA receiver on the APEX telescope. The deep ob- 
servations allow us to detect submm galaxies with observed flux 
densities above ~ 4.5 mJy, while the gravitational magnifica- 
tion reveals galaxies with intrinsic sub-mJy flux densities. We 
derive the magnification of the foreground clusters by using the 
lens equation for clusters modeled as a superposition of Navarro, 
Frenk and White (NFW, 1997) mass density profiles whose pa- 
rameters are inferred from published papers on the selected clus- 
ters. From the magnification maps, we calculate intrinsic flux 
densities and derive submm number counts for the entire survey. 
We carry out a stacking analysis on 24 jum detected sources in the 
fields to probe the correlation between submm and mid-infrared 
emission and detect stacked 870 fim observed flux densities of 
S 870 A<m < 800 yi/Jy for sources that are undetected individually in 
the maps. 

This paper is organized as follows: in Sect.|2]we describe the 
submm observations and data reduction and the Spitzer MIPS 
24 fim archival data; in Sect.[3]we present the resulting maps. In 
Sect. |4] we discuss the lensing models and the number counts 
and in Sect. [5] we present a stacking analysis. Section [6] dis- 
cusses the contribution of our submm signals to the Extragalactic 
Background Light discovered by COBE. The results are summa- 
rized in Sect. |7] 

Throughout the paper, we adopt the following cosmological 
parameters: a Hubble constant Hq = 71 km s _1 Mpc~', a matter 
density parameter Qo = 0.27, and a dark energy density param- 
eter Qao = 0.73. The redshift z = 0.3 where three of our clusters 
reside corresponds to an angular-diameter distance of 911 Mpc 
and a scale of 4.42 kpc/arcsec. z = 2.2, the median redshift of 
known submm galaxies, corresponds to an angular-diameter dis- 
tance of 1728 Mpc and a scale of 8.38 kpc/arcsecQ. 



2. Observations and data reduction 

We have gathered data from galaxy cluster fields observed with 
the L ABOCA bolometer camera on the APEX0 telescope in 
Chile dGiisten et al.ll2.Q06h . The five clusters clusters are merg- 
ing systems, and their high masses yield areas of large gravita- 
tional magnification, which increases the possibility of finding 
very intrinsically dim submillimeter sources lensed by the clus- 
ter. Three of the cluster field observations are from our own ob- 
serving programs, while the AC 114 data (Principal Investigators 



S. Chapman and F. Boone) were downloaded from the ESO 
archive and the Abell 2163 data were provided by the PI. M. 
Nord. Detailed information about the cluster fields, including 
integration time and noise levels of the final maps, is given in 
Table □ 

Ground-based submm observations suffer from the fact that 
the Earth's atmosphere is by far brighter than the astronomi- 
cal sources. The changing temperature of the atmosphere fur- 
ther complicates the data reduction. The relatively small field- 
of-view of LABOCA (11.4') limits the influence of the spatial 
temperature structure of the atmosphere on the measurement. 

LABOCA observes in total power scanning mode, where 
the telescope scans the sky in a pattern that is designed to fa- 
cilitate the retrieval of the astronomical signal and the removal 
of the atmospheric signal. The scanning pattern that was used 
for our observations is an outwards winding spiral which is re- 
peated at four raster points. At a given time during the scan, 
the atmospheric signal is correlated across the entire array and 
we can model and remove it. The faint astronomical signal is 
not correlated, unless it is distributed on scales comparable to 
the fie ld-of-view of the bolometer camera. The Minicrush soft- 
ware (lKovacsll2008ri . that we use to reduce the data, utilizes this 
approach when removing the correlated atmospheric noise. 

Several types of calibration data are taken during the ob- 
servations. Absolute flux calibration is determined from obser- 
vations of the primary calibrators: Neptune, Uranus and Mars. 
When no primary calibrator is available, secondary calibrators, 
which are well studied objects for which the flux ratios to the pri- 
mary calibrators are known, can be used. Measuring the calibra- 
tors also gives a measure of the opacity of the atmosphere. The 
opacity is also measured by performing skydips, which are fast 
scans that measure the sky temperature as a function of elevation 
at constant azimuthal angle. These scans are performed every 2- 
3 hou rs. The calibration o f LABOCA data is described in de- 
tail bv lSiringo et al.l(l2009l) . The telescope pointing was checked 
regularly with scans on nearby bright sources and was found to 
be stable within 3" (rms). The angular resolution (FWHM) of 
LABOCA on APEX is 19.5". 

2.1. Data reduction 



The d ata were reduced using the Minicrush s oftware ([ Kovacs 



l2008h. similarly to the procedure described in IJohansson et al 



1 We used Ned Wright's cosmology calculator ([Wright 200(J avail- 
able at http://www.astro.ucla. edu/~wright/cosmocalc. html. 

2 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory 



d20 1 Oh . We summarize the steps here. The data are organized in 
MBFITS-files, where data from each bolometer as a function of 
time are saved in a so-called timestream. Each scan, and thereby 
each MBFITS-file, contains the timestream data of each bolome- 
ter. Minicrush attempts to remove the correlated noise by tem- 
porarily regarding it as a signal, and fitting a model to all the 
timestreams at the same time. This model is then removed from 
each timestream, and the result is a cleaner signal, with less cor- 
related noise. This procedure is repeated a number of times (for 
LABOCA usually six to eight times) until the resulting signal 
is "white", that is that most of the 1 //-type noise has been re- 
moved. 

An advantage of this method for removing correlated noise 
is that the gains of each individual bolometer can be estimated 
during the process. Another method for determining the gains is 
to observe a bright calibration source and scan it to produce a 
fully sampled map with each bolometer. The information about 
the gains is used to flatfield the data. It can also be used to flag 
and remove suspicious bolometer channels from the reduction. A 
channel with almost no optical response will appear to have very 
low noise level, but searching for and flagging channels with low 
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Table 1. Observed cluster fields. 



Target 


a" [J2000] 


8 a [J2000] 


- 


rms* 


a c 




[h m s] 


[° ' "] 




[mJybeairT ] 


[arcmin 2 ] 


Abell 2163 


16 15 45.1 


-06 08 31 


0.203 


2.2 


150 


Bullet Cluster 1 


06 58 29.2 


-55 56 45 


0.296 


1.2 


220 


Abell 2744 2 


00 14 15.0 


-30 22 60 


0.308 


1.5 


220 


AC 114 3 


22 58 52.3 


-34 46 55 


0.312 


1.2 


130 


MS 1054-03 4 


10 57 00.2 


-03 37 27 


0.823 


1.6 


200 



Notes. <n) Central coordinates of the 2x2 raster square scanning pattern; these positions differ slightly from the central X-ray positions. {b) The 
noise level measured in the central 10 arcminutes of each map, as described in Sect. 13.21 w Extent of the LABOCA maps. 

(1) Alternative name IE 0657-56. Project's observing identification (Obs. ID): O-079.F-9304A-2007, E-380.A-3036A-2007. 

(2) Alternative name AC 118. Obs. ID O-081.F-9319A-2008. 

(3) Alternative name Abell S 1077, Obs. ID E-081.A-0451A-2008, E-078.F-9032A-2007. 

(4) Obs. ID O-083.F-9300A-2009 



gain would find and remove that channel from the reduction. The 
pipeline also flags spikes and glitches in the bolometer channels. 

We used the option ' -deep' in Minicrush. This turns on 
the most aggressive filtering and is useful when searching for 
point-like sources. Extended structures, such as the Sunyaev- 
Zeldovich increment from the clusters, are filtered out. 

2.2. Making maps 

When the pipeline has removed the correlated noise from the at- 
mosphere and from the instrument, flagged any optically dead 
channels and any bad pixels (for example hit by cosmic rays), 
the timestreams should be "white", i.e. free from 1 //-type noise. 
The astronomical signal is typically too weak to be seen in the 
timestreams, and maps from individual scans have to be pro- 
duced and co-added to reduce the noise. The maps are made 
by using the scanning pattern of the telescope and map each 
bolometer position onto a grid of points; when a bolometer has 
"seen" a certain pixel on the map, its flux is deposited there. 
Since several bolometers have seen the same portion of the sky, 
the final flux value in one map pixel is an average of the flux 
of the bolometers that observed that part of the sky, weighted 
by the variance of the individual bolometers. A noisy bolometer 
thus contributes less to the flux density value in a single pixel in 
the map than a less noisy one. 

Together with the flux density map (the "signal" map), a 
noise map is created. The coadded values of the time-stream 
weights are used to create the noise map. A signal-to- noise map, 
which is the signal map divided by the noise map, is also ap- 
pended to the FITS-file. 

2.3. SpitzerMIPS data reduction 

In Sect. |5] we dis cuss a stacking an alysis in the LABOCA maps 
on Spitzer MIPS dRieke et al.1120041) 24 yum source positions. We 
describe here the MIPS data reduction and source extraction pro- 
cedure. We follow the general recommendation not to use the 
pipeline-processed MIPS mosaics, but to re-process the data, 
because persistent images can create artefacts on the mosaics. 
Since we use only the 24 //m MIPS data, "MIPS" will refer to 
that band in the remainder of this paper. 

Four of the clusters in our survey have been observed in the 
24 /mi MIPS band. We used all the available 24 /mi data for each 
cluster field, as summarized in Table|2] We started by download- 
ing the basic calibrated data from the Spitzer science archive. 
For each Astronomical Observation Request (AOR), we created 



a flatfield frame using the script "flatfield_24_ediscs.nl". 
That frame was then used to correct for any persistent problem in 
the data. We then used Mopex to do overlap correction on all the 
data for each target. The overlap-corrected data were then mo- 
saiced using Mopex. We used the default values for all the steps 
in the pipeline. 

Source extraction from the mosaics was performed with the 
Apex tool. We did Point Response Function (PRF) fitting and 
aperture photometry to detect significant MIPS sources. The 
number of detected sources per field is listed in Table _0 to- 
gether with the MIPS coverage across the LABOCA map, the 
median noise level for point sources and the 24 /mi source num- 
ber density. For the Bullet Cluster and MS 1054-03, the MIPS 
maps cover almost the entire LABOCA field, but for AC 114 
and Abell 2744 the MIPS map are significantly smaller. 

Table 2. Summary of archival MIPS data used in this study. 



Cluster name 


£2° 




N s c 






[arcmin 2 ] 


U'Jy] 




[10 3 deg~ 2 ] 


Bullet Cluster 1 


20' X 21' 




325 


2.8 


Abell 2744 2 


6' X 11' 


13.8 


193 


10.5 


AC 1 14 3 


8' X 11' 


10.0 


208 


8.5 


MS 1054-03 4 


17' X 19' 


20.3 


552 


6.2 



Notes. ( '" Angular coverage across the LABOCA field. Median la- 
noise level for extracted 24 jim sources in the part of the map covered 
by LABOCA. (c> Number of sources in the area reported in column 2. 
Program identification numbers (PIDs) 
(1) 40137, 40593 
<2) 83, 3644 
(3 » 83, 50096 

<4) 20740, 83, 3644, 50726 



3. Results 

Figure Q] shows the final reduced signal-to-noise maps of the 
cluster fields. They have been smoothed with a Gaussian of the 
size of the beam (19.5") to a final resolution of 27.5", and emis- 
sion on scales larger than 100" has been filtered out. Contours 
of the noise maps are overlaid on the signal-to-noise maps. 

In the remainder of this section we discuss the source ex- 
traction process, and describe the source catalog and the Monte 
Carlo simulations that we use to characterize the noise level and 
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completeness in each field. We also discuss the method we used 
for flux deboosting and how we estimate the number of spurious 
detections. 



3. 1 . Source extraction 

We impose a detection threshold of > 3.5<x in the signal-to-noise 
maps for source extraction. We exclude sources that lie on the 
edge of the maps where the signal-to-noise representation is not 
accurate. 

In each map we extract any source position with significance 
> 3.5<x and fit a circular, two-dimensional Gaussian to the same 
position in the signal map. We limit the size of the Gaussian 
to that of the beam's FWHM of the images, because submm 
galaxies are ver y likely to be point- like with respect to the 
LABOCA beam. iTacconi et alj (|20Q6.) found a median source- 
sizes of < 0.5", derived from interferometric CO line emission 
observations, in a sample of submm galaxies. In two cases we fit 
elliptical Gaussians, where we know from previous observations 
that the LABOCA sources are comprised of emission from two 
or more galaxies. The two sources are (1) the brightest source 
in the Bullet cluster (SMM J065837.6-555705), which is known 
to b e a blend of two images of the same z - 2.79 galaxy (see 



e.g. iWilson et alj|2008bl: Ijohansson et alj|20ld; iGonzalez et alj 
strongly lensed by the cluster potential and an ellipti- 
cal cluster member; (2) SMM J105656.4-033622 in MS 1054- 
03 which has three SC UBA 850 //m counterparts as reported by 
iKnudsen et all d2008l) . and is discussed further in AppendixlAl 

We measure the noise level for each source from the Monte 
Carlo simulated maps described in the next section. The pro- 
cedure ensures that neither confusion noise nor nearby sources 
contaminate the noise estimate. 

The final source list is displayed in Table[3] There, we list, to- 
gether with positions and measured flux densities of each source, 
deboosted flux densities and gravitational magnification values, 
which are discussed in the following sections. 



3.2. Survey completeness and depth 

Monte Carlo simulations are used to analyze the noise levels of 
the observations and to simulate the completeness of the survey. 
We create so-called jack-knife maps^, which are the result of 
coadding all the scanmaps on one target when multiplying half 
of the scans with -1. This effectively removes any astronomical 
emission from the resulting maps, and makes them representa- 
tions of the instrumental noise only. By randomizing the positive 
and negative scans a large number of different jackknives can 
be created, all of which being random realizations of the noise 
in our observations. For each cluster field we created 500 jack- 
knife maps. 

We note that the confusion noise in the real maps is effec- 
tively removed from the jack-knifed maps. This implies that the 
noise level is underestimated. One can show that at the depth of 
our maps, the instrumental noise exceeds the confusion noise. 
To estimate the confusion level, i.e. the flux level where a larger 
integration time will not decrease the noise level due to the 
unresolved background sources, we use the standard estimate 
that confusion occurs whe n there is one source per 30 beams 
(ICondonl 1 974t lHoggl200 lh . We can estimate the confusion level 
from the relations presented bv lKnudsen et alj d2008l) . They used 



a power law distribution for the number counts (N(> S) = 
NpS with No = 13000 deg~ 2 and a = 2.0) of submm galaxies 
dBarger et al.ll999HBorvs et al.l2003l) . The confusion noise level 

30£2beamM)A< a j where /u is the mean gravi- 
tational magnification across the field. Thus, the confusion level 
is lowered by the lensing, and for the map FWHMs of 27.5" and 
a mean magnification factor of 1 .5 (which is a lower limit), the 
confusion level is < 3 mJy, which is lower than the faintest de- 
tected source in our survey, at 4.6 mJy. The confusion noise is 
thus much smaller than the instrumental noise, and can be safely 
neglected in the following analysis. 



Noise levels 

The jack-knifed maps are used to estimate the noise levels of 
each map as a function of angular distance from the center. In a 
circle of increasing size we extract all pixels in each jack-knife 
map, measure the standard deviation, and then take the average 
of the values for all the jack-knife maps. This procedure is re- 
peated for increasing values of the radial coordinate. The his- 
tograms of the distribution of pixel values in the jack-knife maps 
are well described by Gaussians, whose standard deviations are 
measures of the noise level in each cluster map. In Fig. [2] we 
show the results of that analysis for the five cluster fields. It can 
be seen that the noise level changes very little out to a radius of 
5'. The noise level reached at that distance is the value that we 
report for each map in Table Q] Also, because the noise level is 
almost constant in this region, it is the part of the map that we 
use in the number counts. 

From the jack-knives we also estimate the noise level for 
each detected source. Around the position of the submm source 
(which is not present in the jack-knife maps) we extract a circu- 
lar area the size of the beam, and measure the standard deviation 
for those pixels in each jack-knife map. The average value is 
reported as the uncertainty in the second column in Table [3] 




3 Jack-knifing is a general statistical technique used in all fields of 
science to estimate the precision of sample statistics, and it has been 
used by several groups to analyze mm/submm bolometer data. 



3 4 5 6 
Radius [arcmin] 

Fig. 2. Pixel-to-pixel root-mean-square (rms) as a function of ra- 
dial distance from the map center for the five cluster fields, de- 
scribed in Sect. 13.21 



4 



Johansson et al.: A LABOCA survey of submillimeter galaxies 
Table 3. Flux densities and magnification factors for submm sources detected in the survey. 



Submm source name 


S 87fl iim 
O i\) j-llU. 


S * 

^ deboosted 


" demag 


A* 


S/N'' 




[mJy] 


[mJy] 


[mJy] 






Abeli 2163 












SMM J161525. 8-060803 


7.8 ± 2.2 


4.3 ± 3.5 


3.8 


1.12 


3.5 


*\MM T161S4I 9-060817- 


4 9 + 18 






1.52 


2.8 


SMM TI6I547 7-060948 


8.9 ± 2.0 


6.4 ± 2.3 


3. 1 


2.04 


4.4 


SMM TI6I553 1-060655 


6.7 ±1.8 


4.5 ± 2.3 


2.9 


1.56 


3.8 


IE 0657-56 












SMM J065751.4-560112 


13.5 ±2.7 


9.7 ±3.0 


8.9 


1.09 


5.0 


SMM J0658 13.4-555732 


4.9 ± 1.0 


4.2 ± 1.0 


2.8 


1.54 


5.1 


SMM J065819.4-555830 


8.2 ± 0.9 


7.7 ± 0.9 


5.0 


1.55 


8.7 


SMM J065822.9-560041 


4.8 ± 1.2 


3.8 ± 1.3 


3.0 


1.28 


4.2 


SMM J065824.0-555723 


5.3 ± 0.9 


4.7 ± 1.0 


2.2 


2.10 


5.7 


SMM J065824.5-555512 


15.1 ± 1.0 


14.7 ± 1.0 


9.2 


1.59 


14.9 


SMM J065825.5-555640 


6.9 ±0.9 


6.4 ± 1.0 


2.9 


2.21 


7.3 


SMM J065827. 3-5601 16 


9.0 ± 1.3 


8.0 ± 1.3 


6.4 


1.25 


7.0 


SMM J065828.9-555349 


9.3 ± 1.2 


8.6 ± 1.2 


6.4 


1.35 


7.9 


SMM T06S£3^ 7 SSS441 


4 6 -t- 1 1 
+ .0 :£ 1.1 


6 -t- 1 9 
3.0 ± 1 .Z 


i.j 


1 .JO 


4 1 
+ . 1 


SMM J065837.6-555705 1 


48.6 ± 1.3 


48.0 ± 1.3 


0.6 


75. 3 


36.7 


SMM J065845.6-555848 


6.2 ± 1.1 


5.5 ± 1.1 


3. 1 


1.79 


5.9 


SMM J065846.6-560212 


7.2 ± 1.9 


4.6 ± 2.5 


3.9 


1.18 


3.8 


SMM J065853.2-560046 


7.8 ± 1.5 


6.4 ± 1.6 


5. 1 


1.25 


5.2 


SMM J065853.7-555543 


5.5 ± 1.2 


4.5 ± 1.2 


2.9 


1.55 


4.7 


SMM J065856.0-555652 


5.4 ± 1.2 


4.4 ± 1.3 


3.0 


1.50 


4.6 


SMM J065901. 3-555218 


11.9 ± 2.1 


9.7 ±2.1 


S.5 


1.14 


5.8 


SMM J065915.6-560108 2 


23.6 ±5.9 


— 


— 


1.11 


4.0 


Ahell 2744 












SMM T0014.00 9-309447 


8 1 + 15 


6 8 + 16 


6 2 


1 09 


g 4 


SMM TOO 1406 3-301 Q49 


7.4 +14 


6 1 + 15 


5 


1 23 


^ 2 


SMM J001407.7-302439 


5.8 ± 1.2 


4.8 ± 1.3 


4.1 


1.17 


4.7 


SMM J001408.6-302142 


9.0 ± 1.3 


8.0 ± 1.3 


5.2 


1.55 


6.9 


SMM J001418. 3-3025 25 


4.7 ± 1.3 


3.3 ± 1.6 


2.9 


1.17 


3.6 


SMM J001423.4-302018 


6.9 ± 1.5 


5.5 ± 1.6 


3.5 


1.54 


4.6 


AC 114 












SMM J225835.0-344453 


9.2 ± 1.4 


8.0 ± 1.4 


7.0 


1.15 


6.6 


SMM J225835.7-344812 


5.0 ± 1.5 


3.1 ± 2.0 


2. 1 


1.46 


3.4 


SMM J225844.7-345131 


11.4 ±1.2 


10.6 ± 1.2 


S.3 


1.28 


9.4 


SMM J225905.9-344639 


5.1 ± 1.1 


4.2 ± 1.2 


3.6 


1.19 


4.6 


MS 1054-03 












SMM J105643.7-033543 


6.8 ± 1.5 


5.3 ± 1.7 


5.3 


1.01 


4.4 


SMM J105650.8-034046 


7.9 ± 1.5 


6.5 ± 1.6 


6.4 


1.01 


5.2 


SMM J105656.4-033622 1 


9.8 ± 1.8 


6.5 ± 1.3 


4.4 


1.46 


6.8 


SMM J105703.2-034135 


8.4 ± 1.6 


6.8 ± 1.7 


6.8 


1.01 


5.1 


SMM J105703.7-033309 


8.3 ± 1.7 


6.5 ± 1.8 


6.5 


1.00 


5.0 



Notes. " The flux density is obtained by fitting a beam-shaped Gaussian to the LABOCA source, while the uncertainty is measured in the jack- 
knife maps, as discussed in the text. 

b By construction, the deboosted flux densities have asymmetric uncertainty intervals, but the difference between the upper and lower uncertainties 
is smaller than the number of significant digits reported here. 

c Signal-to-noise ratio of each detection calculated from the observed values in column 2. This value may be smaller than the threshold for source 
extraction (3.5cr) which is imposed in the signal-to-noise map. 

(1) Extended sources. (2) Sources for which the posterior flux distribution from the flux deboosting algorithm has no local maximum at S > 
mJy. SMM J161541. 2-060817 is very close to the detection thresh old, and it is pos s ible th at it is not a real source. It was excluded from the 
number counts calculation. SMM J065915. 6-560108 was discussed in Johansso n et al .1 d201Ct) and has both AzTEC 1.1 mm and Herschel SPIRE 
counterparts. See also the discussion about spurious detections in Sect. 13.41 <3) Our derived m agnification for t his hig hly magnified source is fi ~ 41, 
which is different from the value of 75 which was derived from more detailed modelling by Go nzalez et alj d2009h . and thus we adopt their value. 



Completeness 

We simulate the effects of completeness by inserting artifi- 
cial sources (Gaussians of the size of the LABOCA beam) 
into randomly chosen jack-knife maps, running the source ex- 
traction algorithm, and then comparing the detected sources 
to those input. We hmit the angular area to the central 
10'. Similar analy s is have previously b een performed b y 
Beelen et al.1 (l2008h: iKnudsen et al.l (120081) : IWeifi et all (l2009h : 
Johansson et alj (120 1 Oh . We simulate sources of flux densities 
from 1 to 15 mJy, increasing incrementally by 0.5 mJy, and 
make 500 simulations per flux density bin. Although the jack- 



knife maps are realizations of the noise in the maps, it is possible 
(and consistent with the underlying Gaussian statistics) to find 
fake "sources" which are noise peaks. We therefore include the 
condition that a detected source should be situated sufficiently 
close (within a beam) of the input source. 

The results for the completeness simulations of the five 
LABOCA maps are displayed in Fig. [3] The completeness 
curves follow the general expected behavior; a noisier map has 
a lower completeness value at a certain flux density. From the 
curves we see that for example the Bullet Cluster map is ~ 70% 
complete at 4.2 mJy (the 3.5<x limit for source extraction), while 
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MS 1054-03 at the corresponding flux density of 5.4 mJy the 
map is 65% complete. At a flux density of 6 mJy the maps are 
48% (Abell 2163), 93% (Bullet Cluster), 58% (Abell 2744), 96% 
(AC 1 14) and 77% (MS 1054-03) complete. 

The completeness curves are used to evaluate the submm 
number counts. With the 3.5cr significance limit for source ex- 
traction we must take undetected sources into account when con- 
structing the number counts, which are discussed in Sect. 14.11 
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Fig. 3. Completeness of the five cluster fields as a function of 
flux density, as described in Sect.P 



3.3. Flux deboosting 

The faint population of submm sources acts to boost the 
flux density of the detect ed sources in the s urvey. We have 
used a Bayesian recipe dCoppin et al. 



2005I) to correct for 
this fl ux boosting (see e.g. iScott et al.l I2008 : iHogg & Turner! 
1998b. The procedure is described in detail in Appendix A of 
Johansson et al.l (120 1 Oh . A prior flux distribution is calculated 
in a Monte Carlo simulation, where we create sky maps with 
sources distributed in flux according to a Schechter distribution, 
and source positions are drawn randomly, ignoring the effects 
of clusterings We generate 10 6 simulated maps, and calculate 
a mean flux distribution from the histogram of pixel values for 
each map. This histogram is our prior flux distribution. The prior 
is multiplied with the probability of measuring a flux density S m 
when the intrinsic flux density is S{. This probability is modeled 
as a Gaussian distribution, with the observed flux and noise lev- 
els in column 2 of Table[3]as mean and dispersion. The product 
of the prior flux distribution and the probability to measure a flux 
S m when the intrinsic flux is 5, is normalized to yield the pos- 
terior flux distributiorQ. The deboosted flux density corresponds 
to the .x-axis value found at the local maximum of the poste- 
rior flux distribution. These values are listed in Table [3] Due to 



4 Clustering of faint submm sources have been observed in severa l 
studies (see e.g. lGreve et~aill2004l : rScott et alJl2006l : IWeifi et al.ll2009h . 
but at a resolution of 27.5" the confusion noise contribution from clus- 
tering is much s maller than that from the "normal" Poisson distributed 
confusion noise ( Negrello et al.l|2004t) , and can therefore safely be ne- 
glected. 

3 See the Figure Al in Joha nsson et alj d20 lOh for an illustration of 
these operations. 



the monotonic decrease of the prior flux distribution, the process 
yields non-symmetric flux density uncertainties, but as noted in 
Table [3] the difference between the upper and lower uncertainty 
is smaller than the number of given digits in the table. 

We model the prior distribution by a Schechter function with 
the following form 



£-*(£)--p(-*/n 



(1) 



Schechter! 1976[) w ith parameters from the SHADES survey 



( Coppin et al.1 2006). that have been scaled from 850 pm to 



870 /im using a submm spectral index of 2.7. That yields the fol- 
lowing parameter values: N' = 1703 deg~ 2 mJy -1 , 5" = 3.1 mJy 
and a = -2.0. 



3.4. Spurious detections 

Our adopted criterion for source extraction (S/N > 3.5) means 
that we are detecting sources close to the noise in the maps. It 
is therefore possible that some of our detections are spurious; 
they might be due to a noise peak boosted by confusion noise or 
instrumental artefacts. It is important to investigate how many of 
the sources in our catalog that might be spurious detections. We 
do that by employing two different techniques: 

1. The number of negative 3.5cr peaks in the maps: We run the 
source detection algorithm on inverted maps with the same 
3.5cr criterion, to estimate the number of spurious detections. 
In the five inverted maps we find four sources, indicating that 
at least four of the sources in our catalog could be spurious 
detections. 

2. The probability that a source has a negative deboosted 
flux density: the deboosting algorithm gives us the poste- 
rior flux density for each flux/noise-pair, which is a prob- 
ability distribution for the flux density. For the 37 sources 
in the survey, we find four that have a > 5% probability of 
having a negative flux density, and are possibly spurious. 
These four sources include SMM J161541. 2-0608 17 and 
SMM J0659 15. 6-560 108 which are discussed in the notes 
of Table[3] 

The two methods of estimating the false detections agree well. 
We note that these calculations only gives a statistical measure 
of the number of spurious sources, and not which those sources 
are. However, it is more likely that the least significant sources 
are false detections. 



4. Cluster lens models and number counts 

The clusters in our survey were partly chosen for their lensing 
properties because their high masses lead to areas of high mag- 
nification. In order to estimate the magnification of the detected 
sources, knowledge of the mass distribution of the clusters is 
required. The following calculations use the thin lens approxi- 
mation. 

The magnification factor due to a gravitational lens is given 
by the relation 

u = — — , (2) 
H det^l V ' 

where J\ is the Jacobian of the lens equation, 



m) - do 



1 - K - Jl -y 2 

-Jl 1 ~ K + J\ 



(3) 
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Table 4. Details about the cluster mass models 



Cluster name 


Redshift 


Mass 
110 I4 Mr,l 


Scale radius 
[arcmin] 


c 


Rel. pos. 
[arc sec] 


Reference 


Abell 2163 


0.203 


22 


4.4 


3.6 


0,0 


Radovich et al. (2008) 


Bullet Cluster 


0.296 


31 


4.0 


3.3 


-86, -24 


Clowe et al. (2004) 






8.0 


0.7 


5.3 


86, 24 




Abell I /44 


0.308 


1 1 


2.5 


3.5 


0, 


Boschin et al. (2t)0o) 






3.5 


1.5 


3.9 


48, 135 




AC 114 


0.312 


12 


2.6 


3.5 


0,0 


Camcusano et al. (2001) 






4.3 


1.7 


3.9 


75,-75 






2.3 


1.2 


4.2 


80, 30 




MS 1054-03 


0.823 


3.4 


1.0 


2.9 


0,0 


Hoekstra et al. (2000) 






3.4 


1.0 


2.9 


50, 25 








3.4 


1.0 


2.9 


-60, -20 





k is the convergence of the lens while j\ and j2 are the compo- 
nents of the complex shear. 

We model the clusters as a superposition of one or more 
spherica lly symmetric Nava rro, Frenk & White (NFW) mass 
profiles (iNavarro et al. 1997). For an NFW profile, the conver- 
gence is dTakada & Jainll2003l) 



M m fc 2 

k(6) = r= — = - — 7— — F(c8/8 vir ), 



and the shear is 



r(6) = p^G(ce/e vk ), 



2nr 2 . E c rit 

vir L,lu 

where the critical surface mass density 

D S 



^rrit — 



4ttG D l D ls ' 



(4) 



(5) 



(6) 



where Di is the angular-diameter distance to the lens, D$ the 
distance to the source, and Du the distance between the source 
and the lens, and F, G and / are functions of th e concentration 
param eter which can be found in appendix B in iTakada & Jainl 
(120031) . 

The virial radius of the mass distribution, r vlI can be calcu- 
lated from the virial mass 



3M V 



Pcrit {z)4nA c 



1/3 



where 



Pcrit(z) = H 2 {Z) 

pcrit(O) H 2 



= Q (i +zf + n 



AO 



(7) 



(8) 



is the critical density at the redshift of the lens in a fiat Universe. 
The virial ov erdensity A c can be estim ated from a fit to numerical 
simulations dBrvan & Normanlll998h : 

A c = 187T + 82x- 39x 2 , (9) 

where x = «,„(z) - 1 and a),„(z) = £2o(l + z) 3 H 2 /H 2 (z). 

We estimate the concentration parameter c for a certain mass 
and redshif t using a fit to X-ray luminous clusters of galaxies 
lEttori et alJ ( l2010h 



c(M vil -,z) 



lO 4 / M. 



1 + z \ M» 



(10) 



with M t = 1.0 x 10 15 M o and the fitted parameters A = 0.558 + 
0.008 and B = -0.451 + 0.023. The resulting magnification fac- 
tors and number counts do not depend strongly on the assum ed 
c - M relation; using the relation from lBullock et all (1200 ll) re- 
sul ts in magnification factors that differ from those derived from 
the lEttorietaTld2010t) fit at a level lower than the statistical un- 
certainties. 

We wrote a computer program that generates magnification 
maps by solving Eq. © on a two-dimensional grid. Each clus- 
ter was modeled as one or a sum of NFW halos, whose masses 
were taken from mass models in the literature. The parame- 
ters for these mass models are summarized in Table [4] We do 
not include any individual galaxies in our models. We assume 
the background sources (the source plane) are at a redshift of 
Z = 2.5 but find that the magnifications are not particularly sen- 
sitive to changes in source redshiffl By creating magnification 
maps for the cluster fields, the magnifications of the detected 
sources could be read out from their position in the maps. 

A short discussion of each of the five cluster models follows: 

- Abell 2163: We used the mass from a weak lensing analysis 
performed bv lRadovich et alj d2008l) and modeled the cluster 
as a single NFW profile. Work bv lMaurogordato et aD d2008h 
suggests that the cluster is an ongoing merger and that the 
mass distribution is elongated, which we do not account for 
in our simple model. 

- IE 0657-56 is the most massive cluster in the survey and 
gives the largest area with high magnification of the five clus- 
ters. It consists of two components, o ne main cluster and a 
smaller subcluster. lClowe et al.l d2004l) fitted the main cluster 
to a NFW profile and meas ured the mass of the subcluster us- 
ing ap erture densitometry dClowe et al.ll2000l : iFahlman et al.l 
119941) . 

- Abell 2744 is made up of two sub clusters aligned along the 
line-of-sight dBoschin etalJ 120061) . The subclusters have a 
mass ratio of 3:1 as e stimated from a fit to NFW profiles by 
iBoschin et al.l d2006l) . We use their mass estimate as masses 
for two concentric NFW profiles; 

- AC 114: We use the results from lCampusano et alj d2001l) for 
ou r model. They improve d upon a previous lensing model 
bv iNataraian et alJ dl998l) . Their model of AC 114 is made 
out of a central cluster component, two smaller subclusters 



6 Setting the redshift of the source plane z = 2.0 and z = 3.0 results 
in changes in the magnification factors of less than 10% and on aver- 
age a change of 2%. This uncertainty is smaller than the absolute flux 
calibration uncertainty, or the instrumental noise level for each source. 



7 



Johansson et al.: A LABOCA survey of submillimeter galaxies 



and a galaxy-scale component centered on each bright clus- 
ter galaxy. We include the three large components but not the 
galaxy-scale components in our model. 
MS 1054-03 consists of three distinct mass concentrations. 
We model the cluster as three NFW profiles with masses es- 
timated from a fit of three singular i sothermal sphere profiles 
performed bv lHoekstra et al.l ( 20001) . 



The resulting magnification maps are displayed in Fig. [4] The 
positions of the detected sources for each cluster are overlaid. 
The magnification values are listed in Table [3] as well as the 
demagnified flux densities. In Fig. [5] we show the square root 
of the area (in the image plane) that has a certain magnification 
factor or larger. This figure shows the complex interplay between 
mass and redshift that determines whether there are areas of high 
magnification. Abell 2163, which is the second most massive 
cluster, is a less effective lens than AC 1 14 because it is at lower 
redshift and therefore its mass is distributed over a larger area on 
the sky. 




Magnification n 

Fig. 5. Square-root of the area with magnification > fi as func- 
tion of magnification (//) for the five cluster models in the sam- 
ple. 



4.1. Number counts analysis 

Since gravitational lensing affects both the observed flux den- 
sity of a source and the area surveyed, we have to make cor- 
rections when calculating number counts. The observed flux of 
the sources must be demagnified to estimate the intrinsic flux. 
The intrinsic flux of a source is related to the observed flux by 
Sobs = 1-iSi where is the magnification of the source. Also, be- 
cause the magnification is not constant across the maps, neither 
is the sensitivity. 

We consider the central 10' of the maps where the noise 
is approximately constant (see Fig. [2] and the discussion in 
Sect. 13.21 ). We then impose the same significance criterion as for 
the submm maps: that a source must have a signal-to-noise ratio 
of 3.5 or higher in order to be reliably detected. This signal-to- 
noise ratio corresponds to a certain minimum observed flux den- 
sity S m \ n . A source with intrinsic flux density of 5; is then only 
detected if it lies in a region with magnification /i > S m i n /S i. 

The area of this region in the lens plane, A e ffj(> S), is the 
effective area that we are surveying for sources of a certain in- 



trinsic flux density or greater. The area in the lens plane corre- 
sponds to a smaller area in the source plane due to magnification 
and the effective area we are surveying in the source plane is 



A e ff,s(> s ) = ^A„jU„, 



(11) 



where A„ is the area in the lens plane of a single area element and 
fi„ is the magnification of that particular element. In our case A„ 
corresponds to the area of one pixel and jj n the magnification of 
that pixel. Thus a single detected source corresponds to number 
count of 1 /A e ff tS (> S) sources per unit area. 

We also account for the effects of incompleteness in the 
maps, using the Monte Carlo simulations described in Sect. 13.21 
For each source we detect with a flux density corresponding to 
a completeness of C we expect there to be on average N un i = 
1 /C - 1 undetected sources with the same observed flux density. 

By assuming that those undetected sources are uniformly 
distributed in the map we can calculate the probability that they 
have a certain intrinsic flux. This probability is 



P(S MtlSobs) = 



t/A 



field 



(12) 



where A b s ^int is the area in the image plane which has a magnifi- 
cation in the interval required to place a source with an observed 
flux density 5 b s into the bin corresponding to S mt . Ajj e id is the 
total image plane area of the field in which the source lies. 

4.2. Resulting number counts 
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Fig. 6. Derived number counts for the survey (filled boxes) com- 
pared with results from previous observations. 



The resulting number counts are shown in Fig |6] tog 



fig |pj tog ether 

with the results f rom the surveys of iBeelen et all d2008h and 
IWeifi et all d2009t) carried out with LABOCA at 870 ^m as 
well as those f r om SC UBA surveys of lCoppinet al. ( 2006]) and 
iKnudsen et all (I2008I) at 850 //m, for which the flux values were 
scaled from 850 //m to 870 fim with a spectral index of 2.7. 
There is generally good agreement with previous results. A 
much larger survey of lensing clusters would be needed to reduce 
the uncertainties and make more secure prediction about the 
dim submm galaxies. Uncertainties for the number cou nts were 
calcul ated from Poisson statistics, using the tables in iGehrelsl 
(fT986l) . The number counts and their uncertainties are presented 
in Table 



8 



Johansson et al.: A LABOCA survey of submillimeter galaxies 



Table 5. Number counts. 



S 870 fim 

[mJy] 


N(>S) 
deg - 


N 

1 ' source 


0.6 


9570f"f 

1005 


1 


1.8 


25 5 2 5683 


3 


2.5 


177 1 2797 

1 " 1 1069 


9 


3.5 


7J1 529 


7 


5.0 


501 1334 


2 


6.0 


QT'7567 

-"'188 


7 


7.0 


164 435 


2 



5. Stacking analysis 

We now turn to the analysis of the undetected sources in the 
maps. It is well known that the most numerous contribution to 
the submm galaxy population comes from dim sources with low 
flux densities (see for example the number counts in Fig. |6|. 
Those dim sources cannot be identified individually in the 
submm maps, but their presence can be inferred statistically. 

Stacking (coadding different parts of a map to lower the 
noise) has proven to be an efficient method to reveal the un- 
derlying dim populatio n of submm sourc e s. Such analysis ha s 
been performed by e.g. IScott et all d2008l) : iDunne et alTd2009t) : 
iGreve et alj d20ld) . Except for Abell 2163, all fields in our sam- 
ple have available Spitzer MIPS 24 fim data, as summarized in 
Sect. 1231 

The MIPS catalogs give us positions and flux densities of all 
24 fim sources in the fields, which we use to stack data in our 
LABOCA maps. We excluded MIPS positions which are farther 
than 6' away from the LABOCA map center, where the noise 
level is rapidly increasing (see Fig. [2]). Including positions fur- 
ther out in the map would not lower the noise in the stacked 
signal. We also excluded MIPS positions that lie closer than the 
size of the LABOCA beam from the cataloged submm sources 
in Table [3] For each 24 fim position we extract submaps of size 
2' x 2' from the LABOCA map (labelled St). We also extract 
the same region from the noise map (cr,), and use the following 
relation to stack the submaps: 

Sstack - 2 , (13) 

Z/l/07 

i.e. a summation weighted by the variance. Lastly, we note that, 
although no cataloged submm sources will enter the central 
position of the stacked signal because those positions are dis- 
carded, they may contaminate the outskirts of the stacked map. 
Therefore, we also subtracted models of the cataloged submm 
sources from the LABOCA maps before stacking. This low- 
ers the noise levels in each of the stacked maps but leaves the 
stacked flux densities unchanged. 

In Fig. |7] we show the stacked images for the four clusters. 
Each of the stacked maps shows a significant detection in the 
central region, which is well fitted by a circular Gaussian of the 
size of the beam. Flux densities and noise levels for the maps are 
summarized in Table[6] The noise levels of the stacked maps are 
measured by subtracting the best-fit Gaussians and calculating 
the pixel-to-pixel rms of the residual maps. The flux densities of 
the stacked signals range from ~ 350 - 820 fiiy. This is equal 
to the mean observed flux density of the sources that contribute 
to the stacked signal. In general terms, a deeper MIPS catalog 
yields a lower 870 fim flux value (comparing the stacked fluxes 



Table 6. Results from the stacking analysis. 



Cluster field 


S 870 (im° 

b"Jy] 


[//Jy beam -1 ] 


SJN 


o c 
J demag 




Individual fields 






Bullet Cluster 


815 


88 


9.3 


591 


AC 114 


356 


72 


4.9 


222 


MS 1054-03 


475 


74 


6.4 


413 


Abell 2744 


514 


82 


6.2 


345 




535 


All fields 
37 


14.5 


390 



Notes. (n) Measured by fitting a circular two-dimensional Gaussian 
of the size of the beam to the stacked signal. (6) Pixel-to-pixel rms 
in the residual map where the best-fit Gaussian model was sub- 
tracted. {c) Demagnified stacked signal taking the magnification at each 
24 /jm position into account. 



with the 24 fim depth the deeper maps have a lower stacked sig- 
nal). See also the following section. 

The stacked signal corresponds to an observed 870 fim flux 
density. To investigate the intrinsic fluxes of the dim galaxies we 
use the magnification maps derived from the cluster models. We 
find the magnification factor for each 24 fim position in the map 
and then stack the submaps extracted from the LABOCA map 
again, this time dividing each submap by the magnification of 
the central source. This calculation is only valid for the central 
part of the stacked map, since the magnification is not constant 
across the submaps. We fit a circular Gaussian to the stacked 
signals, and report the measurements in the fifth column of Table 
[6] The demagnified stacked fluxes are lower than the original. 

Finally, we coadd the four stacked signals, weighted by their 
noise-maps, to find the total stacked 24 fim signal for the entire 
survey (excluding Abell 2163). The map, shown in the fifth panel 
in Fig. [7] is a 14.5<x-detection, as reported in Table [6] When 
performing the same operation on the stacked maps corrected 
for gravitational magnification, we find a mean signal of 390 fiiy 
for the four cluster fields. What would be the properties of an 
submm galaxy wi th such a flux density ? Assuming a median 
redshift of z — 2.2 dChapman et alj2005l) . dust temperature T t j ~ 
40 K and dust emissivity index a — 2, the 870 fim flux density 
of 390 fiJy corres ponds to a far-infrared l uminosity Lfir ~ 6.4 x 
10 11 L (eq. 11 in lDe Breuck et al.ll2003h . By assuming that the 
submm emission originates mainly from starburst phenomena, 
which follows a Salpeter initial mass function with a low-mass 
cutoff mi = 1 .6 M Q , we can estimate a sta r -form ation rate SFR ~ 
6OM yr~ 1 using Eq. 4 of lOmont et al.1 d200ll) . It is clear that 
the stacking analysis uncovers a population of submm sources 
different from that detected directly in the maps. The derived 
far-infrared luminosity and star-formation rate depend on several 
assumptions about the underlying submm population. 

Several groups have obtained si milar results when stacking 
on MIPS or radio source positions. IScott et all d2008l) found a 
stacked signal of 324 + 25 fiiy on ~ 20 00 MIPS positions i n 
the AzTEC study of the COSMOS field dScoville et alj|2007b . 
Their MIPS map had a similar depth to those in the present 
study. Accounting for the wavelength difference by scaling the 
AzTEC 1 . 1 mm flux density with a mm/submm spectral index 
of 2 - 3, the 870 fim flux density would be ~ 520 - 660 fiiy. 
This is s lightly higher than the intrinsic LABOCA flux density, 
390 fiJy. IGreve et al.1 (EoToh found similar stacked flux values in 
their study of the Extended Chandra Deep Field South. 



9 



Johansson et al.: A LABOCA survey of submillimeter galaxies 



5.1. On the validity of the stacking detections 

In order to assess the validity of the detection, we perform the 
stacking analysis on random positions in the LABOCA map, 
drawn from a uniform distribution within the same map area 
as the 24 /an maps. We run 20 such simulations for each field. 
None of the simulated maps has stacked signal with a signifi- 
cant 870 pm source in the center. The nondetection in the sim- 
ulated maps gives confidence in the stacking results on the real 
24 pm position, and shows a correlation between the MIPS and 
LABOCA maps. We discuss the nature of this correlation in 
more detail in the next section. 
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Fig. 8. Top panel: normalized histogram of the signal-to-noise 
values in the LABOCA maps at the 24 pm positions (dash-dotted 
line), compared with the mean histogram from 20 sets of random 
positions (solid line). Lower panel: Difference signal between 
the signal-to-noise histogram above, showing that the main part 
of the stacked signal is due to sources with signal-to-noise ra- 
tios between 0.5 and 2. This histogram clearly shows a deficit of 
points at negative signal-to-noise units and an excess at positive 
signal-to-noise units. 



To investigate the significance of the sources that contribute 
to the stacked signals, we extract pixel values at the 24 pm posi- 
tions from the signal-to-noise maps and compare them to the ran- 
domly distributed positions. The histograms for all the 24 pm po- 
sitions not ascociated with a significant LABOCA source is 
shown in Fig. [8] It is compared with the histogram for ran- 
domly distributed points, which has the shape of a Gaussian. 
The difference signal between the two curves indicates that 
LABOCA points with significance 0.5cr < S/N < 2cr contribute 
the most to the stacked signal. 



5.2. The 524/an -Sgvo/™ relation 

So far, we have investigated the signal resulting by stacking 
LABOCA sub-map at each MIPS position within the fields. 
However, it is not plausible that all MIPS sources contribute 
equally to the stacked map. By choosing subsets of the total 
MIPS catalog with different 24 pm flux values, we examine a 
possible correlation between the flux density of the stacked sig- 
nal and the 24 pm flux density. We perform this analysis in the 
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Fig. 9. Results from the stacking analysis in MS 1054-03, 
AC 114 and Abell 2744, as described in Sect. EH Boxes: mea- 
sured flux density in the maps when dividing the 24 pm po- 
sitions into six equal parts, and stacking them. The stacked 
flux densities are measured by fitting a circular Gaussian to 
the maps, except in the lowest two 24 pm flux bins where no 
significant stacked signal was detected, and we instead mea- 
sured the flux density in a circular aperture with the diameter 
of the beam FWHM. Circle: the stacked signal for all the three 
clusters. Shaded bars: histograms indicating the distribution of 
24 ^m flux densities within each of the six flux bins. Each of the 
bins have two bars and the different shades of gray discriminate 
between them. In the highest flux bin the rightmost bar includes 
all the 24 pm sources with flux densities larger than 650 pJy. All 
flux values were corrected for gravitational magnification. 



fields of MS 1054-03, AC 114 and Abell 2744. There are 918 
24 pm sources in the three fields. We divide the catalog into six 
sub-catalogs with equal number of sources (138), with median 
24 pm flux densities of 72, 121, 164, 207, 303 and 599 pJy, and 
perform the stacking analysis for each of them. Both the 24 and 
870 pm fluxes are demagnified. We find significant signals in the 
four highest 24 pm flux bins. The results are plotted in Fig. [9] At 
low 24 pm flux (S24fim< 300 pJy) we find a linear relation be- 
tween the 24 pm and the 870 pm flux. At higher flux densities a 
turnover occurs and the c urve flattens out. T he results are along 
the same lines as those of iGreve et al.l d2010l) who find a flatten- 
ing of the S24/jm - S&j0fj m relation at 524^ m ~ 350 piy, in their 

sta cking analys is. 

IGreve et all d2010t) argue that the linear relation at low MIPS 
fluxes is an indication that those sources are dominated by 
star formation, whereas the flattening of the curve at larger 
24 pm fluxes is due to contamination by active galactic nuclei 
(AGN). The mid-IR flux is sensitive to warm dust, which is 
likely to be heated by an AGN. The 870 pm flux is more sen- 
sitive to colder dust, heated by starbursts. While the mid-IR flux 
increases the 870 pm flux stays constant, because it is not sensi- 
tive to the warm dust emission. 

5.3. A MIPS source contributing to the stacked signal 

It is beyond the scope of this paper to present a full multi- 
wavelength analysis and comparison between the MIPS and 
LABOCA maps. Here, we note one source which has a large 
magnification and is almost detected in the LABOCA map. 
iRigbv etai] ( l2008l) presented Spitzer/IRS spectroscopy of lensed 
galaxies, and discussed one source in the center of AC 114, 
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gravitationally magnified by a factor of 9.7, and at a redshift of 
Z = 1 .47. At this position in our LABOCA map there is a positive 
signal with a significance of 3.4<x and flux density ~ 4 mJy. 

Figure[lO]shows a postage-stamp cutout of the region around 
the source. The LABOCA counterpart to the Spitzer source is 
just below our detection limit, and we would need more data 
to confirm it. If the 870 /mi source is real, the high magnfication 
value would mean that its intrinsic flux density S &jo IJm < 0.5 mJy. 
We note that our lensing model for AC 1 14 gives a magnification 
value of yu = 6 for the Spitzer source position. Given that it lies 
very close to the brightest cluster galaxy in AC 1 14, the small 
discrepancy between the two magnification values are likely due 
to the lack of modeling of individual cluster member galaxies in 
this work. 
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Fig. 10. Postage stamp cutout of the core in AC 1 14. The circle 
(with diameter of the angular resolution of the LABOCA map, 
FWHM = 27.5") deno tes the posit i on of the Spitzer source at 
Z = 1 .47, as reported bv lRigbv etai](l2008h . 



6. Contribution to the extragalactic background 
light 

We have described and quantified the submm emission from 
two types of objects in the paper: the significantly detected 
sources in the maps, and those whose mean flux density is in- 
ferred from the stacking analysis. We now estimate the frac- 
tion of the Extragalactic Background Light (EBL) detected in 
the far-infra red by COBE dPuget et al.lll996tlFixsen et al.lll998t 
lDweketal.iri998b that our survey has resolved. At 870 jum the 
surfac e brightness of the EBL is 44 + 15 Jy deg 2 dGreve et alJ 
l20loh . 

The total flux density of the detected submm sources that 
lie within the 10 arcminute central region is ~ 250 mJy. This 
corresponds to a surface brightness of 2.3 Jy deg -2 , thus showing 
that the LABOCA observations have resolved ~ 5% of the EBL 
into significant > 3.5<x-sources. Another way to calculate the 
contribution to the EBL is to use the number counts, which have 
been corrected for completeness. Integrating the number counts 
yields a surface brightness of 5.7 Jy deg -2 , corresponding to ~ 
13% of the EBL. 

Turning to the stacked signal, where we detected a mean 
observed flux density of 535 + 37 fiJy at 1278 positions in an 
area of ~ 500 arcmin 2 , the surface brightness is ~ 4.9 Jy deg -2 
and the stacked submm signal thus corresponds to 1 1 % of the 
EBL. Thus, in total our observations have uncovered the source 



of 24% of the EBL. Because gravitational lensing preserves sur- 
face brightness we can choose to perform this calculation using 
either the observed or intrinsic flux densities and areas. 

Other au t hors find different EBL contributions. 
iKnudsen et al.l d2008l) find that their observations resolve 
almost all of the extragalactic background light, since they 
discovered seven galaxies with sub-mJy i ntrinsic flux levels , 
and thus probe the number counts very deep. lGreve et alJ (f20loh 
find that the contributions for the stacked signal in the Extended 
Chandra Deep Field-South varies with redshift between 10% at 
z = 0.5 to 40% at z = 2. 

7. Conclusions 

We used the LABOCA receiver on APEX to carry out a submm 
survey of five clusters of galaxies. The clusters act as gravita- 
tional lenses and magnify background sources. The main results 
of the survey are summarized below. 

1. We discovered 37 submm sources, out of which 14 are new 
submm detections. 

2. We modeled the galaxy clusters as superpositions of spher- 
ical NFW halos and generated magnification maps for the 
five clusters. 

3. The magnification maps were used to correct for the gravita- 
tional lensing and to obtain the intrinsic flux densities of the 
detected sources. 

4. We constructed number counts taking into account both the 
gravitational lensing and the varying completeness level. The 
number counts probe the sub-mJy level and are consistent 
with previous work within the uncertainties. 

5. We performed a stacking analysis in the LABOCA maps on 
positions of detected 24 //m sources in the fields. The stack- 
ing yields > 4.9cr detections in all fields with MIPS cover- 
age, reaching noise levels below 100 //Jy, more than an order 
of magnitude deeper than the individual maps. 

6. By dividing the 24 /urn catalog in MS 1054-03, AC 1 14 and 
Abell 2744 into six equal halves, we find a linear relation 
between S 24 fan and 5870 at low 24 fim fluxes, followed by 
a flattening of the relation at 524^ m ~ 300 /vJy. This behavior 
can be explained if the low MIPS fluxes trace star formation 
while the higher values are dominated by AGN heating. 

7. The observations reveal a total of ~24% of the infrared extra- 
galactic background light, where ~13% comes from the sig- 
nificant submm sources and ~11% comes from the stacked 
signal. 

Acknowledgements. We would like to thank the APEX staff for excellent sup- 
port during the observations. We also thank Rodrigo Parra for his help during 
observations of Abell 2744 in June 2008. We thank Martin Sommer for supply- 
ing us with his data on Abell 2163 and Eiichi Egami for letting us use his cata- 
log of MIPS source in the Bullet Cluster field. We would like to thank Kirsten 
Knudsen for helpful discussions. We acknowledge support from the Swedish 
research council (Vetenskapsradet). We thank the referee for very useful com- 
ments. 

This work is based in part on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with NASA. This research has made use 
of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. This research has also made 
use of NASA's Astrophysics Data System Bibliographic Services. 

References 

Austermann, J. E., Dunlop, J. S., Perera, T. A., et al. 2010, MNRAS, 401, 160 
Barger, A. J., Cowie, L. L., & Sanders, D. B. 1999, ApJ, 518, L5 



11 



Johansson et al.: A LABOCA survey of submillimeter galaxies 



Beelen, A., Omont, A., Bavouzet, N, et al. 2008, A&A, 485, 645 
Bertoldi, R, Carilli, C, Aravena, M., et al. 2007, ApJS, 172, 132 
Blain, A. W., Smail, I., Ivison, R. J., Kneib, J.-R, & Frayer, D. T. 2002, 

Phys. Rep., 369, 111 
Borys, C, Chapman, S., Halpem, M., & Scott, D. 2003, MNRAS, 344, 385 
Boschin, W., Girardi, M., Spolaor, M., & Barrena, R. 2006, A&A, 449, 461 
Bryan, G. L. & Norman, M. L. 1998, ApJ, 495, 80 
Bullock, J. S., Kolatt, T. S., Sigad, Y., et al. 2001, MNRAS, 321, 559 
Campusano, L. E., Pello, R., Kneib, J., et al. 2001, A&A, 378, 394 
Chapman, S. C. Blain, A. W., Smail, I., & Ivison, R. J. 2005, ApJ, 622, 772 
Chapman, S. C, Scott, D., Borys, C, & Fahlman, G. G. 2002a, MNRAS, 330, 

92 

Chapman, S. C, Smail, I., Ivison, R. J., & Blain, A. W. 2002b, MNRAS, 335, 
L17 

Clowe, D., Gonzalez, A., & Markevitch, M. 2004, ApJ, 604, 596 

Clowe, D., Luppino, G. A., Kaiser, N., & Gioia, I. M. 2000, ApJ, 539, 540 

Condon, J. J. 1974, ApJ, 188, 279 

Coppin, K., Chapin, E. L., Mortier, A. M. J., et al. 2006, MNRAS, 372, 1621 
Coppin, K, Halpem, M., Scott, D., Borys, C, & Chapman, S. 2005, MNRAS, 
357, 1022 

Danielson, A. L. R., Swinbank, A. M., Smail, I., et al. 2010, MNRAS, 1565 

De Breuck, C, Neri, R., Morganti, R., et al. 2003, A&A, 401, 91 1 

Dole, H., Lagache, G., Puget, J., et al. 2006, A&A, 451, 417 

Dunne, L., Ivison, R. J., Maddox, S., et al. 2009, MNRAS, 394, 3 

Dwek, E., Arendt, R. G., Hauser, M. G., et al. 1998, ApJ, 508, 106 

Egami, E., Rex, M., Rawle, T. D., et al. 2010, A&A, 518, L12+ 

Ettori, S., Gastaldello, F., Leccardi, A., et al. 2010, A&A, 524, A68+ 

Fahlman, G., Kaiser, N, Squires, G., & Woods, D. 1994, ApJ, 437, 56 

Fixsen, D. J., Dwek, E., Mather, J. C, Bennett, C. L., & Shafer, R. A. 1998, ApJ, 

508, 123 
Gehrels, N. 1986, ApJ, 303, 336 

Gonzalez, A. H., Clowe, D., Bradac, M., et al. 2009, ApJ, 691, 525 
Gonzalez, A. H., Papovich, C, Bradac, M., & Jones, C. 2010, ApJ, 720, 245 
Greve, T. R., Ivison, R. J., Bertoldi, E, et al. 2004, MNRAS, 354, 779 
Greve, T. R., Wei/3, A., Walter, E, et al. 2010, ApJ, 719, 483 
Glisten, R., Nyman, L. A., Schilke, P., et al. 2006, A&A, 454, L13 
Hoekstra, H., Franx, M., & Kuijken, K. 2000, ApJ, 532, 88 
Hogg, D. W. 2001, AJ, 121, 1207 
Hogg, D. W. & Turner, E. L. 1998, PASP, 110, 727 
Ikarashi, S., Kohno, K., Aguirre, J. E., et al. 2010, ArXiv e-prints 
Johansson, D., Horellou, C, Sommer, M. W., et al. 2010, A&A, 514, A77+ 
Knudsen, K. K., Barnard, V. E., van der Werf, P. P., et al. 2006, MNRAS, 368, 
487 

Knudsen, K. K, van der Werf, P., Franx, M„ et al. 2005, ApJ, 632, L9 
Knudsen, K. K, van der Werf, P. P., & Kneib, J.-P. 2008, MNRAS, 384, 1611 
Kovacs, A. 2008, in Presented at the Society of Photo-Optical Instrumentation 

Engineers (SPIE) Conference, Vol. 7020, Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series 
Maurogordato, S., Cappi, A., Ferrari, C, et al. 2008, A&A, 481, 593 
Natarajan, P., Kneib, J., Smail, I., & Ellis, R. S. 1998, ApJ, 499, 600 
Navarro, J. F, Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 
Negrello, M., Hopwood, R., De Zotti, G., et al. 2010, Science, 330, 800 
Negrello, M., Magliocchetti, M., Moscardini, L., et al. 2004, MNRAS, 352, 493 
Nord, M., Basu, K., Pacaud, E, et al. 2009, A&A, 506, 623 
Omont, A., Cox, P., Bertoldi, E, et al. 2001, A&A, 374, 371 
Patanchon, G., Ade, P. A. R., Bock, J. J., et al. 2009, ApJ, 707, 1750 
Perez-Gonzalez, P. G, Egami, E., Rex, M., et al. 2010, A&A, 518, L15+ 
Puget, J.-L., Abergel, A., Bernard, J.-P, et al. 1996, A&A, 308, L5+ 
Radovich, M., Puddu, E., Romano, A., Grado, A., & Getman, F. 2008, A&A, 

487, 55 

Rex, M., Ade, P. A. R., Aretxaga, I., et al. 2009, ApJ, 703, 348 
Rex, M., Rawle, T. D., Egami, E., et al. 2010, A&A, 518, L13+ 
Rieke, G. H., Young, E. T., Engelbracht, C. W., et al. 2004, ApJS, 154, 25 
Rigby, J. R., Marcillac, D., Egami, E., et al. 2008, ApJ, 675, 262 
Schechter, P. 1976, ApJ, 203, 297 

Scott, K. S., Austermann, J. E., Perera, T. A., et al. 2008, MNRAS, 385, 2225 
Scott, K. S., Yun, M. S., Wilson, G. W., et al. 2010, MNRAS, 405, 2260 
Scott, S. E., Dunlop, J. S., & Serjeant, S. 2006, MNRAS, 370, 1057 
Scoville, N, Aussel, H, Brusa, M., et al. 2007, ApJS, 172, 1 
Siringo, G., Kreysa, E., Kovacs, A., et al. 2009, A&A, 497, 945 
Smail, I., Ivison, R. J., & Blain, A. W. 1997, ApJ, 490, L5+ 
Smail, I., Ivison, R. J., Blain, A. W., & Kneib, J.-P. 2002, MNRAS, 331, 495 
Swinbank, A. M., Smail, I., Longmore, S., et al. 2010, Nature, 464, 733 
Tacconi, L. J., Neri, R., Chapman, S. C, et al. 2006, ApJ, 640, 228 
Takada, M. & Jain, B. 2003, MNRAS, 344, 857 

Vieira, J. D., Crawford, T. M., Switzer, E. R., et al. 2010, ApJ, 719, 763 
Wardlow, J. L., Smail, I., Wilson, G. W., et al. 2010, MNRAS, 401, 2299 
WeiB, A., Kovacs, A., Coppin, K., et al. 2009, ApJ, 707, 1201 



Wilson, G. W, Austermann, J. E., Perera, T. A., et al. 2008a, MNRAS, 386, 807 
Wilson, G. W., Hughes, D. H, Aretxaga, I., et al. 2008b, MNRAS, 390, 1061 
Wright, E. L. 2006, PASP, 118, 171 1 

Zemcov, M., Borys, C, Halpern, M., Mauskopf, P., & Scott, D. 2007, MNRAS, 
376, 1073 



Appendix A: Comparison with other submm 
observations 

MS 1054-03 

Even though the galaxy clusters in our sample have been well 
studied across the electromagnetic spectrum, sub mm maps have 
been p ublished for only one system, MS 1054-03 dKnudsen et al.l 
120081 hereafter K08) using SCUBA. S ubsets of those SCUBA 
data were previously analyzed by Chapman et all l2002at 
IZemcov et alj [2007; Kn udsen et al J 12005 , with mostly similar 
source catalogs. Differences in the SCUBA source catalogs 
are discussed by K08. We will now compare our catalog for 
MS 1054-03 with theirs. 

K08 reach a noise level of 0.86mJybeam _1 in the deep- 
est part of the map, while their area-weighted noise level is 
1.49mJybeairf '. They detect nine significant sources in their 
map. The noise level in the LABOCA map in the central part of 
the map is 1.3 mJy beam -1 (see Fig. [2]), while the average noise 
level within the central 10' is 1 .6 mJy beam 1 . The SCUBA map 
covers the cluster region (14.4 arcmin 2 ) while the usable map 
area in the LABOCA map is ~ 150 arcmin 2 . 

Only one of our detected sources lies in the area covered by 
K08. It is extended with respect to the LABOCA beam, and has 
an angular size of 30" x 35". Its flux density is 9.8 +1.8 mJy. 
In the same area, K08 report three sources, separated by of 
25.3", 25.8"and 18.0". LABOCA, with a coarser resolution 
than SCUBA, causes the three sources to blend together. We 
compared the measured LABOCA source size with a simple 
model of the three K08-sources, constructed as a sum of three 
LABOCA beam shaped Gaussians (angular FWHM 19.5"). We 
then smoothed this model image with a Gaussian of angular 
FWHM of the size of the LABOCA beam, and then fitted an 
elliptical Gaussian to the resulting map, similarily to what was 
done in the real LABOCA maps. The fitted source has an angular 
size of 33" x 40". The slightly larger angular size of the model 
compared to the observed LABOCA source can be explained 
with uncertainties in the fitted FWHMs and in the SCUBA and 
LABOCA positions. The sum of the flux density of the three 
SCUBA sources, scaled from 850 fim to 870 fim with a submm 
spectral index of 2.7, is 11.9 + 1.5 mJy. This is within the lcr 
uncertainty interval of the LABOCA flux measurement. 

The other six sources detected by K08 have 850 fim flux den- 
sities that, when extrapolated to 870 fim, are too faint to be de- 
tected in the LABOCA map. 



Abell 2163 



iNord et alJ d2009h presented the first LABOCA map of a 
galaxy cluster detected in the Sunyaev-Zeldovich increment, 
Abell 2163. A bright point source close to the cluster center, with 
a flux density of 1 1 .9 + 1.9 mJy, was n oted, but not discussed. 
Using the same data set as in lNord et all but filtering out most of 
the extended SZ-signal, we detect the same point source in our 
map, with a flux density of 8.9 + 2 mJy. 
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1 E 0657-56 

The brightest submm so urce in the Bullet Clu ster has been 
thoroughly discussed (see iJohansson et al.ll2010l and references 
therein). Submm and FIR observations by the SPIRE and 
PACS instruments on the Herschel satellite show counterparts of 
LABOCA sources dRex et al.ll20ToHPerez-Gonzalez et all2010l) 
as part of the Herschel Lensing Survey dEgami et al .II2010I) . 
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Fig. 1. Signal-to-noise maps of the five cluster fields. White circles represent the significant sources in the map and black contours 
show the noise maps of each cluster map at levels of 2, 4 and 8 mJy beam . The signal-to-noise representation causes the appearance 
of the increasing noise towards the edge of each map to be suppressed. 
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Fig. 4. Magnification maps for the five clusters. The positions of the detected submm sources are marked on the maps with circles. 
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Fig. 7. Stacked 870 /im maps (in units of mJy beam - ) on the 24 fim positions for the cluster fields with MIPS observations, overlaid 
with signal-to-noise contours. The white contours range between 3cr and 13cr with an increment of 2<x, while the black countours 
show the -3<t level. The maps have not been corrected for gravitational magnification. The fifth map is the coadded signal of the 
four individual stacked maps, which yields a 14.5<x-detection. In Table[6]we present fits to the stacked maps. 
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